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We engineered a two-dimensional magnetic lattice in an elongated strip geom¬ 
etry, with effective per-plaquette flux ^4/3 times the flux quanta. We imaged 
the localized edge and hulk states of atomic Bose-Einstein condensates in this 
strip, with single lattice-site resolution along the narrow direction. Further, we 
observed both the skipping orbits of excited atoms traveling down our system’s 
edges, analogues to edge magnetoplasmons in 2-D electron systems {1,2), and 
a dynamical Hall effect for bulk excitations (3). Our lattice’s long direction 
consisted of the sites of an optical lattice and its narrow direction consisted of 
the internal atomic spin states (4,5). Our technique has minimal heating, a fea¬ 
ture that will be important for spectroscopic measurements of the Hofstadter 
butterfly (4,6) and realizations of Laughlin’s charge pump (7). 
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Introduction In materials the quantum Hall effects represent an extreme quantum limit, where 
a system’s behavior defies any description with classical physics. In the modem parlance, the 
integer quantum Hall effect (IQHE) for two-dimensional (2-D) electronic systems in magnetic 
fields (8) was the first topological insulator (9): a bulk insulator with dispersing edge states — 
always present in finite-sized topological systems— which give rise to the IQHE’s signature 
quantized Hall resistance (10). 

In classical systems the magnetic field acts purely through the Eorentz force, while in quan¬ 
tum systems a particle with charge g in a uniform field B additionally acquires an Aharonov- 
Bohm phase 0 AB/ 27 r = AB after its path encircles an area A normal to B. (Here <l>o = 
2nh/q is the flux quanta and 2nh is Planck’s constant.) We engineered a 2-D lattice with square 
geometry for ultracold charge-neutral atoms in which we directly controlled the acquired phases 
as atoms traversed the lattice, giving a tunneling phase 0AB/27r ^4/3 around each plaquette. 
These phases take the place of the Aharonov-Bohm phases produced by true magnetic fields and 
suffice to fully define the effective magnetic field. Aharonov-Bohm phases of order unity — 
only possible in engineered materials {11,12), or in atomic (5, 13-17) and optical {18) settings— 
fragment the low-field Eandau levels into the fractal energy bands of the Hofstadter butterfly {6). 
Such Hofstadter bands are generally associated with a non-zero topological index: the Chern 
number {10). 

Topologically nontrivial bulk properties are reflected by the presence of edge channels, com¬ 
posed of edge states, with quantized conductance. In fermionic systems, the number of edge 
channels is fixed by the aggregate topological index of the filled bands (9, 10,19)', this ulti¬ 
mately gives rise to phenomena such as the IQHE for electrons. Conceptually the constituent 
edge states can be viewed as skipping orbits (9, 20, 21): in the presence of a strong magnetic 
field, nascent cyclotron orbits near the boundary reflect from the hard wall before completing 
an orbit, leading to skipping trajectories following the system’s boundary. In contrast, localized 
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bulk states correspond to closed cyclotron orbits. 


By applying large effective fields to atomic Bose-Einstein condensates (BECs), we directly 
imaged individual, deterministically prepared, bulk and edge eigenstates. In IQHE systems 
these states would govern the conductivity, but as individual eigenstates they exhibit no time- 
dependance. The corresponding dynamical entities are edge magnetoplasmons, consisting of 
superpositions of edge eigenstates in different Eandau levels (i, 2), or here magnetic bands. We 
launched these excitations and recorded their full motion for the first time, observing both a 
chiral drift along the system’s edge and the underlying skipping motion. 

The Harper-Hofstadter Hamiltonian 



( 1 ) 


with complex hopping amplitudes tx and tg governs the motion of charged particles moving 
in a 2-D lattice (<5, 22) with sites labeled by j and m, the situation which we engineered for 
our neutral atoms. Analogous to the Eandau gauge in continuum systems, we describe our 
experiment with real tg (no phase), and with complex tx = \tx \ exp(—dependent on m. 
As shown in Eig. [IP^ the sum of the tunneling phases around any individual plaquette is 0 ab- 
We implemented a 2-D lattice by combining a conventional optical lattice to define the 
long axis of our system (Oa; direction, with sites labeled by j), with three sequentially coupled 
internal “spin” states to define the short axis of our system (e^ direction, with just three sites 
labeled by m G {—1,0, +1}); in parallel to the work described here, an analogous scheme has 
been realized for fermionic ytterbium (5). As in Eig.[Tp, this system effectively has an infinite 
repulsive potential for \m\ > 2, allowing for the formation of robust edge states. In each band 
of our engineered lattice (Eig.[Tp;), the momentum along specifies the position in e^, (denoted 
by color on the curves), just as for 2-D electrons in Eandau levels. 
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Technique As shown in Fig.[T]4-B, we used BECs in the / = 1 ground state hyperfine 
manifold (23), eonfined in an optieal dipole potential from two 1064 nm laser beams aligned 
along Oa; and e^,, with trap frequeneies (cja;, Uy, uJz)/27: = (50,40,110) Hz. We adiabatieall)|^ 
loaded sueh BECs onto the sites of the 1-D optieal lattiee formed by a pair of = 1064.46 nm 
laser beams eounter-propagating along e^,. This laser defines the single-photon reeoil momen¬ 
tum hki = 27 :%/Xl and reeoil energy = %^k\/2m-^^a-, where m^b is the atomie mass. The 
lattiee depth V ^ 6El gave a hopping strength |ta,| 0.05El along e^,. 

We then eoupled the three \m) sub-levels either with two-photon Raman transitions or an 
rf magnetie field. The Xr = 790.04 nm Raman lasers also eounter-propagatec^^ along e^,, with 
wave-veetor kR = ±27 i/Xr; the rf eoupling effeetively had kR = 0. Either field gave a labora¬ 
tory tunable effeetive tunneling strength \ts\ ~ |fj.| along the spin “direetion” e^, proportional 
to the Rabi frequeney Hr. The 2%kR momentum imparted from these transitions resulted in 
a spatially periodie phase faetor exp{i2kRx) = exp(z0ABj) aeeompanying the ehange in m, 
where 0AB/2vr = kR/ki ±4/3 for Raman eoupling and 0AB/27r = 0 for rf eoupling (4,24). 
The sign of c/xb was eontrolled by the relative detuning of the Raman lasers (see Eig. [T]4-B). 

Eigure shows the band strueture for our system featuring the three Hofstadter bands ex- 
peeted for 0AB/27r ~ 4/3, with our boundary eonditions. The system is most easily understood 
after making the loeal gauge transformation |j, m) —>■ exp (—i^ABjrn) |j, m) whieh transfers 
the Peierls phase from tg into This fully maps our system to Eq. (1), and as shown in Eig. [IP 
the result is an effeetive flux <h/<ho = 0AB/27r per plaquette. 

Experiment We began our experiments by direetly imaging adiabatieally loaded eigenstates 
of the ground Hofstadter band with either 0AB/27r = 0 or 4/3, in an isotropie lattiee with 

'We linearly increased the lattice potential from 0 to ~ in 200 ms. 

^The Raman lasers counter-propagated at an angle 9 = 0.75° away from e^, effectively reducing by « 1 x 
10“^ and introducing a quadratic Zeeman like term with strength —0.0017?^. Both of these effects are negligible 
in our experiment. 
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\tx\ ~ \ts\- After preparation, we used a measurement proeedure eommon to all experiments: 
we simultaneously removed all potentials and eoupling fields (tos < 1 /ws), which returned the 
atoms to bare spin and momentum states. The atomic cloud expanded for a 18 ms time-of- 
flight (TOP) period. During TOP a 2 ms magnetic gradient pulse was applied, Stern-Gerlach 
separating the three \m) states. The resulting 2-D column density was recorded using standard 
absorption imaging techniques, giving the normalized momentum distributions nm{kx) with 
perfect single lattice site resolution along e^. 

Pigure contains typical data with 0 ab = 0, where we adiabatically loaded into the 
ground state and observed nm{kx). The fractional population Um = f (^kxnm{ka) resembles 
(Pig. 1^) that of a particle in a discretized box along e^, while the momentum distributions, 
typical for atoms in an optical lattice (25), have the same profile for each m site. This demon¬ 
strates that the two directions are uncoupled at 0 ab = 0. 

The data in Pig. |^-E, with (/)ab 7^ 0, are qualitatively different as a function of both k^ 
and m. These differences can be understood in analogy with a 2-D electron in a perpendicular 
magnetic field, confined in one dimension with hard walls. Along the confined direction the 
wavefunction is localized to the scale of the magnetic length = \Jk/qB, with center position 
at kxi% in the bulk, where hkx is the electron’s canonical momentum. Por large \kx\, the elec¬ 
tron becomes localized near the edges, lifting the degeneracy of the otherwise macroscopically 
degenerate Landau levels. Each of these points finds an analogue in our observations (Eig. 

Eig. I^-H). In our system the magnetic length = \j2)/2n 0.7 in units of lattice period is 

of order unitjj^ significantly narrowing the bulk state (Eig. [^) as compared to the 0 ab = 0 
case (Eig.[^). In addition, we observed the appearance of states localized at the system’s edges 
(Eig.[3^,H), which are completely absent when 0 ab = 0. These localized edge states are the 
analogue to the current carrying edge states in fermionic IQHE systems. 

^Since all physics is 27r periodic in the acquired phase, our flux 0AB/27r « 4/3 is equivalent to c/ab ~ 27r/3. 
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Having described the static properties of this system, we now turn to dynamics. We loaded 
our systems on to different m sites with ts = 0, and then abruptl)|^ turned on tg, allowing tun¬ 
neling along e^. The resulting initial states all consisted of coherent superpositions of magnetic 
band eigenstates with crystal momentum which began to coherently tun¬ 

nel along and experienced an associated Lorentz force along e^.. Atoms initialized on the 
bulk m = 0 site demonstrated a dynamic Hall effect. Those starting on the edge sites became 
cold-atom analogues to edge magnetoplasmons: they began cyclotron orbits, were reflected 
from the hard wall, and skipped down one edge or the other. 

The dynamics of atoms initialized in the bulk (on the m = 0 site) are presented in Fig.|^ As 
schematically illustrated in Fig. and plotted in Fig. [^, a balanced population oscillated in 
and out of the the originally empty m = ±1 sites as a function of time r. When 0ab 7^ 0, this 
motion drove transverse, i.e. Hall, edge currents Jm=±i(r) = nra{T){Vm) along (Fig.[^), 
where ( 14 „) is the mean velocity of atoms on site m along e^. 

As shown in Fig.|^, a chiral current X = Ji —/_i developed with overall sign following that 
of 0ab- As atoms tunneled to the edges they acquired a transverse velocity controlled by two 
parameters: 0ab set the crystal momentum acquired while tunneling, and gave the natural 
unit of velocity 2tx/hkL. This led to the observed in-phase oscillation of X and the combined 
m = ±1 populations (|m|). 

This synchronous oscillation implies a linear dependance of X on (|m|) whose slope we 
label S, plotted by the red and black symbols in Fig.[^. We confirmed the system’s chirality by 
inverting 0 ab and verifying that S changed sign. For comparison, we repeated the experiment 
with 0AB = 0, and observed no chiral current (empty symbols). These data are all in good 

For the function of estimating the magnetic length - a continuum concept - it is suitable to use (.*g « 

"^We commanded linear ramps to the Raman lasers giving a 300 /rs tum-on time. This time scale was adiabatic 
with respect to the « 5El band spacing in the 1-D optical lattice, yet nearly instantaneous with respect to the 
« Atx magnetic band width. Still, this non-zero turn-on time gave rise to the small offset in Fig. Up, which is 
correctly predicted by our full theory. 
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agreement with our theory (eurves) using parameters obtained from fits to Fig.|^. 

The dependence of the chiral current on the tunneling anisotropy tg/tx shown in Fig.[^ is 
reminiscent of the optical lattice experiments in Ref. (2(5). The chiral current remained linear 
in (|m|), with slope S essentially constanj^in units of 2tx/hkL', this demonstrates a new ki nd 
of dynamic Hall effect. In contrast, the peak edge current X^ax (pink dashed line in Fig. |^) 
strongly depends on tg/tx (Fig. |^), increasing from zero then reaching saturation. For small 
tg/tx, few atoms tunneled giving a correspondingly small X ^a x, while as tg/tx increased Xmax 
began to saturate as essentially all atoms tunnelled (27). 

We then moved from studying bulk excitations to edge excitations by launching edge mag- 
netoplasmons: superpositions of edge states across magnetic bands with crystal momentum 
Qx/ki = T^ab/tt. We created them on either edge, with the potential tilted along as in 
Fig.|^-B, such that the initially occupied site was at the potential minimum. In Fig.|^-D, we 
plot the time evolving average position (m(r)) along and velocity (vx{t)) = Y^m along 
ex- Data shown in pink/blue solid circles are for initial sites (m(r = 0)) = ±1, both of which 
evolved periodically in time but with opposite velocities. The spatial trajectories are illustrated 
in Fig. 1^, where we obtained the displacement {6j (r) ) by directly integrating {vx (r) )/a, where 
a = Al/ 2 is the lattice period. These data clearly show edge magnetoplasmons with their chiral 
longitudinal motion, and constitute the first experimental observation of their edge localization 
and transverse skipping motion. 

Conclusion We realized a large artificial magnetic field — with 0AB/27r 4/3 flux per 

plaquette — in a 2-D lattice geometry. The use of internal atomic states as one of the lattice di¬ 
rections made single site “spatial” resolution along that axis straightforward, enabling our direct 

^The slight downward curvature in our theoretical model, plotted by the grey curve, evidences 5%-level cor¬ 
rections proportional to jig ^ to the tight binding model [Eq. (1)], owing to mixing with excited Bloch bands of the 
1-D optical lattice. 
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observation of chiral edge states, a dynamical quantum Hall effect, and edge magnetoplasmons. 
This and related approaches (5) have the technical advantage over other techniques for creating 
artificial gauge fields in that minimal Raman laser coupling is required (typically 10 to 50 times 
less than previous experiments using Raman coupling (28)) thereby minimizing heating from 
spontaneous emission and enabling many-body experiments which require negligible heating 
rates. Lifetimes from spontaneous emission with this technique are in excess of 10 seconds, 
while all other approaches for creating large artificial gauge fields have lifetimes well below 1 
second (14-16,29). 

With our hard-wall potential, a realization of the Laughlin charge pump (7) is straightfor¬ 
ward: as particles accelerate along e^;, mass moves from one edge to the other in the orthogonal 
direction e^. Remarkably, extending our technique to periodic boundary conditions —coupling 
together the \m = ±1) states— should produce systems exhibiting a fractal Hofstadter spec¬ 
trum (4), even given a three-site extent along e^. Going beyond conventional condensed matter 
realities, the flexibility afforded by directly laser-engineering the hopping enables the creation 
of Mobius strip geometries: topological systems with just one edge (30). 
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Figure 1: Hybrid 2-D lattice. (A) N ce 10^ atom ®^Rb BECs, confined in a 1-D optical lattice 
with a bias magnetic field B^e^, were illuminated either by a pair of counter-propagating Raman 
lasers or an rf magnetic field B^i. (B) The Raman or rf field linking the three internal \m) states 
with Rabi frequency Hr was detuned from the gfisBo/h ^ 0.817 MHz or 1.35 MHz Zeeman 
splitting by 5. The corresponding quadratic Zeeman shift additionally lowered \m = 0) by 
e = 0.05El or 0.1377^. The Raman lasers’ relative phase was actively stabilized at a beam- 
combiner adjacent to the optical lattice retroreflection mirror. (C) The lattice along can be 
considered as a square well with hard walls at m = ±2 for which d ^ 0 tilts the potential 
towards one wall or the other. (D) The 2-D hybrid lattice, where the non-spatial dimension is 
built from the internal states \m) with an effective magnetic flux per plaquette <h/<l)o = 4>ab /27r. 
(E) Three lowest magnetic bands (rainbow colors), computed for our full lattice without making 
the tight binding approximation, with parameters V, 5, e) = (0.2, 6,0, 0.05)£'i. The pale 

curves were computed for RHr = 0. 
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Figure 2: Adiabatically loaded eigenstates. Left panels: Site-resolved, normalized momentum 
distributions nm{kx) obtained by absorption imaging. Right panels: fraetional population rim- 
(A-B) Atoms loaded into the single ground state present for 0 ab = 0 show the expeeted sepa¬ 
rable behavior along e^, and e^. (C-H) Atoms loaded into the upper edge, bulk, and lower edge 
states present for 0AB/2vr 4/3 demonstrate the eoupling between and and loealization 
along e^. 
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Figure 3: (A) Development of ehiral edge eurrents. A system prepared at site m = 0 at time 
Ti = 0 obtains a ehiral edge eurrent at T 2 whieh returns to zero at r^. (B) Fractional population 
versus time for atoms initialized at m = 0 at time r = 0. The dynamics for i^AB are the 
same. (C) Chiral edge current X versus time. Data shown in red, black, and empty circles are 
taken for positive, negative, and zero 0ab, respectively. (D) X plotted against (|m|). The solid 
curves (theory) use parameters V, 5, e) = (0.73, 6, 0.001, 0.05)ii^A determined from B 

for 0 ab 7^ 0, giving tg = O.IAEl, and {hVlji, V, 5, e) = (0.47, 6, —0.01, 0.05)El for 0 ab = 0. 
(E) Maximum edge current versus asymmetry (tg/tx). Inset: slope S (taken from data as in D) 
is nearly independent of tg/tx. 
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Figure 4: Skipping orbits. (A-B) Schematic of dynamics starting from m = — 1 including a 
nonzero detuning which tilted the lattice along e^. (C) Mean displacement (m) versus time r 
for excitations on both edges. Data for systems initialized on m = 1 (m = —1) are depicted by 
pink (blue) circles. (D) Group velocity along e^. versus time for both edge excitations. Numer¬ 
ical simulations (solid curves) use parameters V, 5, e) = (0.58, 5.2, ±0.087, 0.13)£'i 

from fits to population dynamics like C. In both cases, the potential gradient along from 5 
was selected so the initial site (m = ±1) had the lowest energy. (E) Edge magnetoplasmon 
trajectories, where the displacement {Sj{T)) was obtained by integrating {vx{t))/ a. 
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